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Abstract 

We re-investigate the evolution of the strongly degenerate neutrinos in the early universe. With 
the larger degeneracy, the neutrino number freezes at higher temperatures because the neutrino 
annihilation rate decreases. We consider very large degeneracy so large that the neutrino number 
freezes before events in which the particle degrees of freedom in the universe decrease (e.g. the 
muon annihilation and the quark-hadron phase transition). In such a case, the degeneracy by the 
time of nucleosynthesis becomes smaller than the initial degeneracy. We calculate how much it 
decreases from the initial value on the basis of the conservation of the neutrino number and the 
total entropy. We found a large drop in the degeneracy but it is not large enough to affect the 
current constraints on the neutrino degeneracy from BBN and CMBR. 
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I. INTRODUCTION 



It is well known that the evolution of the universe deeply depends on the properties 
of neutrinos (see Ref. for review). In the standard cosmology, we assume three types 
of massless neutrinos with the same number of particles and antiparticles. In this paper, 
we deviate from the last assumption and consider cosmological effects caused by neutrino- 
antineutrino number asymmetry. We use a terminology "degeneracy" for this asymmetry. 
We also assume, just for simplicity of the description, the neutrino number is larger than the 
antineutrino number (i.e. the neutrino has positive chemical potential). When considering 
the opposite case, we need to just exchange the role of neutrinos and antineutrinos. 

Since degeneracy increases the sum of the energy density of neutrinos and antineutri- 
nos, it significantly affects standard predictions such as the big bang nucleosynthesis and 
the cosmic microwave background (in addition, electron-type neutrinos destroy neutrons 
so the nucleosynthesis depends strongly on their number density). Naively, these effects 
monotonously increase with the degree of degeneracy so the observations can put upper 
bounds on the degeneracy in each generation of neutrinos. 

However, things are not as simple as this because the neutrino number freeze-out tem- 
perature also increases with the degeneracy (the rise in the freeze-out temperature is caused 
by the scarcity of the antineutrinos, the annihilation partners of the neutrinos). The compli- 
cation occurs when the freeze-out takes place before the muon-antimuon annihilation ends. 
In this case, when muons and antimuons annihilate, already frozen neutrino number can 
not change but the entropy is transfered to the neutrinos through the elastic scattering 
with electrons and positrons so the neutrinos keep the same temperature with the rest of 
the cosmic plasma. In order to conserve the number while the temperature changing faster 
than the inverse of the scale factor, the neutrino degeneracy parameter, chemical potential 
divided by temperature, decreases. As a result, by the onset of nucleosynthesis, the neutrino 
energy density becomes much lower than the value calculated using the initial degeneracy. 
Therefore, there is a logical possibility that very large initial degeneracy cannot be ruled out 
by the observations. 

This possibility is pursued in Ref. || but with an incorrect picture on the "neutrino 
decoupling" when very large degeneracy exists (the same lines of argument are still found 
in a few papers [||] |7j ||). They have regarded the neutrino number freeze-out and its 



2 



kinematical decoupling from the rest of the cosmic plasma occur simultaneously. The correct 
picture is the one as described in the previous paragraph, the kinetic equilibrium holds well 
after the chemical equilibrium ceases to hold. The reason is that there is small number of 
annihilation partners, antineutrinos, but elastic scattering partners, electrons and positrons, 
are abundant. This is pointed out in Ref. but they have not calculated the relation 
between the initial degeneracy and the final degeneracy after the entropy producing events. 
Obtaining this relation is the main purpose of this paper. 

In the next section, we calculate the neutrino number freeze-out temperature and show it 
increases exponentially with the initial degeneracy parameter. We also justify the neutrinos 
and anti-neutrinos are in kinetic equilibrium during the annihilation by demonstrating their 



elastic scattering rate is sufficiently large. In Sec. |T|, for a certain range of initial degener- 
acy parameter, we calculate their final value after the entropy producing events using the 
neutrino number conservation and the total entropy conservation. This is our main result. 



In Sec. IV, we summarize the results and discuss the current constraints on the neutrino 



degeneracy are not affected. 

II. NEUTRINO NUMBER FREEZE-OUT 

We calculate the neutrino number freeze-out temperature Tf by 1^(7/) = H(Tf), where 
Ta is the rate of the change in the neutrino number density n v through the neutrino- 
antineutrino annihilation processes and H is the cosmic expansion rate. 

In order to calculate T^, we have to sum all the annihilation process rate. Assuming 
one type of the neutrino v\ has degeneracy, we consider the annihilation to some fermion- 
antifermion pairs, Viijpi) +^(^2) — > F{Pz) + F{p4) where the variables in the brackets denote 
the four-momentum of each particle in the comoving frame. For F, there are electrons and 
two non-degenerate types of neutrinos. In addition, we consider the annihilation to muon- 
anitimuon pairs when T > and to quark- ant iquark pairs when T > Tqcd where we 

assume the quark-hadron phase transition to occur instantaneously at Tq CD = 200 MeV and 
the quark-gluon phase contains u, d and s quarks. The contribution of this process to the 
r A is, 
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1 f d 3 p 1 d 3 p 2 d 3 p 3 d 3 p 4 - 2 

x(2n) 4 5^(p 1+P2 -p 3 -p 4 )f Ul (E 1 )f l7l (E 2 )[l - f F (E 3 )][l - f F {E A )}. (1) 

F's are well-approximated to be massless so the square of the invariant matrix element 
\M.{vi + vi — > F + F)\ 2 can be written in the form 32G F [b(pi-p 3 ) 2 + c(prp 4 ) 2 ] where 
G F = (292.80GeV)~ 2 is the Fermi coupling constant. For F ^ I, only the neutral cur- 
rent contributes so that b = (Cy — C^) 2 and c = (Cy + C^) 2 and for F — I, the charged 
current contributes in addition so that b = (C l v — C l A ) 2 and c = (C l v + C^ + 2) 2 . The vector 
and axial- vector couplings (Cf , Cj) are (1/2,1/2) for F — u, (-1/2, -1/2 + 2 sin 2 0) for 
F = e and /i, (1/2, 1/2 - (4/3) sin 2 6) for F — u, and (-1/2, -1/2 - (2/3) sin 2 0) for F = d 
and s where the weak-mixing angle sin 2 9 = 0.231. /'s are the distribution functions of the 
particle species on the subscript and we use the equilibrium form, f FF = l/[exp(E/T) + 1] 
and f UuUl = l/[exp((E±fi)/T) + 1] where fi is the chemical potential and the signs are — for 
ui and + for z/ ; . Here, we assume /i = \i n = — /v P expecting the chemical equilibrium holds 
initially at very high temperature. Some of the integration can be performed analytically 
and we obtain 

q2 f-oo f-oo t-1 rf + 

n^Uu^FF = t8 -^e J o dx J o d y J dz J dt fvi( x )M(y) x3 y 3 (i -z) 2 s~ 5 

x \b + c) {35 4 + 3(x - y) 2 (2t - x - y) 2 - (x - y) 2 5 2 - (2* - x - y) 2 5 2 } 
_ 4(6 - c )(x- y)(2t -x- y)5 2 \ [1 - f F (t)][l - f F (x + y-t)] 
= T 8 G 2 F L(0 (2) 

where x = E 1 /T, y = E 2 /T, z = cos9,t = E 3 /T, S(x, y, z) = (x 2 + y 2 + 2xyz) 1 ' 2 and /± = 
(x+y±5)/2. With these variables, the distribution functions are /^(x) = l/[exp(x±£) + l] 
and f FyF (x) = l/[exp(x) + 1], where we define the degeneracy parameter £ = fi/T. Finally, 
the degenerate neutrino number density is 

^(0 = ^/^T^^M(0T3. (3) 

The cosmic expansion H is determined by the total energy of the universe by the Einstein 
equation, 

[S^ / jggg T 2 
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where M P = 1.2 x 10 19 GeV is the Planck energy and we have introduced g* to express 
the total energy density as p to t = fi , *(O( 7I " 2 /30)T 4 . On calculating the total energy, we do 
not have to worry about temperature dependence of g* because it turns out that when 
the freeze-out temperature becomes so high as muons to contribute as relativistic degree of 
freedom, £ should be large making the muon (and the other particle species other than the 
degenerate neutrino) energy density negligible compared to that of the degenerate neutrino. 
The same is true for pions and the quark-gluon plasma. Therefore we include the contri- 
bution from photons, electrons, positrons, two types of non-degenerate (anti)neutrinos and 
the degenerate (anti)neutrinos, 



0.(0 = 9 
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(5) 



e x ~t + 1 e x '+€ + 1 
This gives usual value g*(0) = 43/4. 

Combining Eqs. ([!]) ~ (H) gives freeze-out temperatures of the neutrino number as a 
function of £ 

"M(0 lV3 



(6) 



7)(£) = 0.999756 g^f' 6 Vj^\ MeV. 

The results are shown in Fig. |I[ We find that Tf increases exponentially with £. This 
is expected because the number density of the antineutrinos, the annihilation partners of 
neutrinos, is exponentially suppressed due to the degeneracy. To see this more explicitly 
and for the later convenience, we here derive the approximate expression for Tf when £ 
is large. First, Eq. (g) is simplified by setting fuXu) ~ exp(— y) exp(— £) and neglecting e ± 
Pauli blocking factors. Using the integration approximation formula for a function satisfying 
<p(0) = o, 

(p'(x)d,X , ^ 7T 2 



we obtain 



ip' (x) f v {x)dx = 

HO 
M(0 



o exp(x - + 1 



¥?(o + ^"(o + S^ {4) (o + 
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For <?*(£), we neglect the exponentially suppressed contribution from antineutrinos and ob- 
tain 

4 4 V 71 J V 7T 



(10) 
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FIG. 1: The relation between the neutrino degeneracy and its number freeze-out temperatures. 
The solid, dotted and small-dotted lines respectively show when v e , v„ and v r has the degeneracy. 
They basically increase exponentially with the degeneracy. At T ~ and Tqqd, according to 

the variation in the relativistic freedom of the fermions, muons and quarks, to which neutrinos can 
annihilate, the change in the freeze-out temperature with degeneracy is retarded. 

We have the expression for Tf by putting these approximations into Eq. (§) and it reproduces 
closely the numerical result plotted in Fig. [I] when £ > 5. 

Now we explicitly write down the expression for Tf taking the leading order in £ to 
compare with that of Refs. [[[]], @ and (4|. Our analysis shows 

T f « 1.37408 (J^\ Vexp(£/3)MeV. (11) 

This and theirs all agree with respect to Tf oc exp(£/3) which originates from the exponen- 
tially suppressed degenerate n v . Our expression contains the factor £ 1 ^ 3 , which originates 
from the degenerate v distribution affecting both of n v and H, but it seems not to appear 
in the others. However, the expression of Ref. |l| agrees with ours if thermally averaged 
momentum divided by T, (p/T), is calculated using the degenerate v distribution so that 
(p/T) « (3/4)£ (they use (p/T) « 3 which is true when the degeneracy is small) and #*(£) 
is approximated as oc £ 4 . As for Ref. [|J, Tf oc £ -2 ^ 3 but this is thought to originate from the 
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incorrect division by n e when deriving Ta- This does not make difference when calculating 
Tf with no degeneracy but is not appropriate when strong degeneracy exists. Ref. has 
not found such a factor of power of £ but this is thought to originate from calculating the 
annihilation rate in the center-of-mass frame. In addition, they also calculate ri e /n e so, as 
Ref. has pointed out, their result has to be corrected anyway. 

Until now, we have not cared about the possibility of the neutrino spectral distortion and 
have used thermal distributions to calculate various quantities. We now justify this proce- 
dure. Since neutrinos with larger energies have more probability to annihilate (from Eq. (^j), 
we see f v / f u oc Ei), they seems to freeze-out later and the neutrino spectrum would distort 
very much. But as we discuss next, the neutrino-electron elastic scattering occur sufficiently 
so the thermal distribution of neutrinos is preserved at the time of annihilation. Since the 
presence of neutrino large degeneracy ensures that the neutrino-antineutrino annihilation 
process occurs with their distribution preserving equilibrium form, it is considered that the 
estimation of neutrino number freeze-out temperature from Ta = H turns out to be good. 

We investigate how these degeneracy parameters undergo a change in the next section. 
Before we proceed, we estimate the temperature at which the neutrinos decouple kinemati- 
cally from e . This temperature is expected to be much lower than the number freeze-out 
temperature because the partners of the elastic processes, electrons and positrons are abun- 
dant regardless of the neutrino degeneracy. So they are coupled until relatively low temper- 
ature and have the same temperature with the others even after their number has frozen out 
at higher temperature. This decoupling temperature is well estimated from the rate for the 
elastic scattering processes because the energy transfer between e and v occurs efficiently 
in this case. The energy transfer efficiency is estimated as follows. The center-of-mass frame 
differencial cross section for the elastic processes is 

(P^) +(t^) = T fG^ + Cl){4 + (l + cos^}, (12) 

\dcos9 J „ +e __ w+e _ \dcos9J v+e+ ^ +e+ 16tt 

where s and 9 are the total initial energy squared and the polar angle in the center-of-mass 
frame. The energy transfer from to v with initial energy (in the laboratory frame) E e ± 
and E v is 

-^transfer = \{ E e± ~ E v){ 1 ~ COS0). (13) 



From (12) and (|Tq) , the expectation value of the energy transfer by one collision is computed 



7 




FIG. 2: The decoupling temperatures of the degenerate neutrinos. The neutrinos are in kinetic 
equilibrium with and T v follows T e ±'s variation above these temperatures. Note that they 
increase with £ but do not exceed 
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Therefore, roughly speaking, v becomes almost as energetic as after one collision. 
Now, we compute the rate for the elastic scattering processes 

* dn ") 1/n. (15) 

and search a temperature Td at which it becomes equal to the expansion rate, just as we 
have done in the case of the annihilation process. The result is shown in Fig. |]. We see 
that Td increase with £ but rather slowly. Therefore, as is the case of no degeneracy, the 
neutrinos are in kinetic equilibrium holding the same temperature as the rest of the plasma 
until at least the muon-antimuon annihilation ends. 
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III. THE CHANGE IN CHEMICAL POTENTIAL AND ENERGY DENSITY 



In this section, we study how the neutrino or its degeneracy parameter £ evolves after 
its number has frozen out. £ is conserved when the annihilation occurs frequently so that 
chemical equilibrium holds. Even after the annihilation practically ceases, £ does not vary 
as long as the temperature T falls as T oc a" 1 where a is the cosmo logical scale factor. But 
this condition is not satisfied when, for example, the muon annihilation progresses because 
the temperature falls more slowly (recall from the end of the last section that the neutrino 
temperature evolves along with the e ± temperature because they are in kinetic equilibrium). 
For the case when the freeze-out takes place before the muon annihilation ends, there is a 
period in which the neutrino number at the freeze-out temperature is conserved while the 
temperature falls slower than a -1 . This causes £ to decrease during the muon annihilation 
or other entropy producing processes. 

We can compute the value of £ after the muon annihilation on the basis of neutrino 
number and total entropy conservation after the freeze-out. We briefly explain why the 
latter holds. The second law of thermodynamics for the total gas in physical volume V is 



where p is the energy density, p is the pressure and n is the number density. The first two 
terms on the right hand side vanishes according to the total energy conservation. For the last 
term, there could be contributions from particle species with non-zero chemical potential, 
but d(nV) = for neutrinos because it has been already frozen out and d(nV) is negligible 
for antineutrinos because its number density is exponentially suppressed (rip = (T 3 /ir 2 )e^ 
for £ — 1) for such large degeneracy as to make the neutrino number freeze out before the 
muon-antimuon annihilation ends. 

The total entropy conservation is expressed as 



where the subscript "/" denotes the value at the neutrino number freeze-out and "final" 
after the p^ annihilation but before the e ± annihilation. £ is conserved when the chemical 
equilibrium holds so we have written initial instead of £/. s to t is the sum of entropy density 
s = (p+p — pn)/T of all the particle species exist. On the right hand side, there exists only 



1 



TdS = d(pV) + pdV - pd(nV), 



(16) 



a 





(17) 
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the relativistic particles, photons, e ± 1 and neutrinos one type with degeneracy and 2 types 
without, so the T 3 can be scaled out as 

s t JT,Q « |> (g,(0 - / ^Y-b) - (18) 

where <?*(£) is the effective degrees of freedom same as the one appeared in Eq. (]5|). We 
neglect exponentially suppressed antineutrino contribution. When we calculate the left hand 
side of Eq. (|i~T|), finite masses of muons (m^ = 106 MeV) and pions (m n = 135 MeV) are 
included in order to treat the annihilation of these particles continuously. We assume the 
quark-hadron phase transition to occur instantaneously at Tqcd = 200 MeV and the quark- 
gluon phase contains u, d and s quarks which are well approximated as massless. Then 



\ 2tt 2 

+4/ M + (T) + 3/-(T)e(T QC D - T)J = — J(T,0T 3 , 



(19) 



where 9{x) is 1 for x > and is otherwise. If denotes the contribution from the massive 
particles i as 



4tt 4 



00 x 2 I x 2 

dx -. — r \ x 2 + af H ^=^^_ 

exp (y^ 2 + a 2 ) ± 1 V 3v/x 2 + a 2 



(20) 



where a, = uii/T and the sign is +/— when z is a fermion/boson. 
The neutrino number conservation is 

a^n v (Tf (initial), ^initial) = a final n f (^filial, £final) , (21) 

and n u is calculated by Eq. (^). Dividing each sides of Eq. ( |T7| ) by Eq. (^TJ), we obtain an 
equation with respect to £fi na i for given initial, 

J(Tf (f initial), 6nitial)M(e fi nal) ~ #(&nal)M(£ initiaJ ) = 0. (22) 

We numerically solve this for £ nna i in the range < initial < 15. The result is shown in Fig. || 
and we see considerable difference between the initial £„ and the final £ w . 

We translate this results in terms of neutrino energy density as shown in Fig. |] (with 
regard to only the degenerate family). We show the final value of the energy density of the 
degenerate neutrinos plus antineutrinos divided by the one with no degeneracy (so-called 
effective number of neutrino types minus two). 
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FIG. 3: The relation between the degeneracy which exists initially (initial) an d which remains after 
the muon-antimuon annihilation (£final)- For initial ^ 6, since the neutrino number freezes while 
the universe contains relativistic degrees of freedom of the particles other than photons, e and 
neutrinos, ^fi na i becomes smaller than initial- 

For a comparison, in Fig. |5|, we also show the degenerate neutrino final energy densities 
when some assumptions are different. We compute the case when the degeneracy parameter 
£ maintains the initial value and the case when we assume, as in Refs. || |§ , the neutrinos 
are not heated by annihilation processes after its number has frozen out. Our result demon- 
strates that for some range of the degeneracy, the final neutrino energy density can be lower 
than the one with the initial degeneracy but not as low as previous analysis [|3| has shown. 

At last, we make a comment on the evolution of the antineutrino degeneracy parameter, 
£ p . At the initial stage, the chemical equilibrium holds so £p = This relation could 

be broken after the neutrino number freeze-out because the neutrino is not in the chemical 
equilibrium any more. However, since the frequent elastic scattering with electrons ensures 
neutrinos and antineutrinos to be in kinematical equilibrium and annihilation is still efficient 
for antineutrinos, the relation £ p = —£ v continues to hold. Therefore, as £„ decreases during 
for example the muon-antimuon annihilation, £ p follows its variation with the sign opposite. 
Meanwhile, when the temperature drops down to the electron-positron annihilation, the 
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FIG. 4: The final energy density of the degenerate type of neutrinos and antineutrinos normalized 
to the one with no degeneracy: p u+ p(^)/ p u+ p(0). 

elastic scattering is not efficient to keep v (and also v if the degeneracy is not very large) 
in the kinetic equilibrium with e . In this case, the momentum distribution is distorted 
from the equilibrium form so, strictly speaking, the notion of the chemical potential or 
the degeneracy parameter is lost. But the distortion can usefully expressed as momentum 
dependent chemical potential and in this sense, + £ p ^ as shown in Ref. where the 
evolution of the neutrino spectrum is fully simulated with £ < 1. 

IV. CONCLUSION 

We have re-investigated some properties of the thermal history of the early universe with 
very large neutrino degeneracy. We have justified and adopted the correct picture of the 
neutrino number freeze-out as explained by Ref. PJ. We have made some complementary 
arguments and calculations to those found in Ref. [Q] and have obtained the results con- 
cerning the evolution of the strong neutrino degeneracy and in turn its energy density. We 
find there are cases that the neutrino degeneracy parameter £ becomes smaller than the 
initial value after the muon annihilation. However, they do not seem to require alteration 



12 



800 




Sinitial 



FIG. 5: A comparison with the results which appear under different assumptions (for u e ). The solid 
line is our result. The dotted line is the energy density calculated from the initial degeneracy. The 
dot-dashed line is the one expected when the neutrinos are assumed to be kinematically decoupled 
simultaneously with its number freeze-out. 



to the cosmological bounds on the neutrino degeneracy such as —0.01 ^ £ Ve ^ 0.22 and 
l&v-l ^ 2-6 || (more stringent bound, < 0.07 for all three neutrino types, is likely 
to apply taking into account neutrino oscillations with maximal mixing HlQ"| . The analysis 
on how the bounds are modified for the region of the Large Mixing Angle solution to solar 
neutrino problems is found in Refs. JL1] and [f[2]j), because, in Fig. [|, the local minimum of 
the final degeneracy caused by the quark-hadron phase transition does not have the value 
lower than the present upper bound. 

Finally, we note our results stem from the fact that the degenerate neutrinos are in kinetic 
equilibrium with the rest of the cosmic plasma well after their number has frozen out. The 
same is true for other possibly degenerate species of stable particles other than the neutrinos 
and our analysis is applicable to them. So, although our analysis on the neutrinos has turned 
out not to affect the cosmological bound on their degeneracy, the remark made in this paper 
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is worth bearing in mind. 
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